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Abstract 



I \ We analyze the influence of different on-site and inter-site interactions on 

the Curie temperature of transition metal magnetic elements. The numerical 

calculations show the well known result that the on-site Coulomb repulsion 

^2 ■ helps ferromagnetism mainly at half-filling of the band, where the 3d elements 

are antiferromagnetic (Cr, Mn). The inter-site interactions, which in the 

'^ \ new approximation decrease the width of the band, favor ferromagnetism 

^ \ at both ends of the band. At the same time they lower the ferromagnetic 

Q ' Curie temperature towards experimental values, removing the paradox which 

persisted for a long time. 
T-H ! 

^ : PACS codes: 75.10.Lp; 75.50.Bb; 75.50. Cc 
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1. Introduction 

Recently there has been a resurgence of research [1-8] into the influence of on-site and 
inter-site interactions on band magnetism. The authors of this paper have shown how 
different inter-site interactions would affect the ferromagnetism in the presence of the on- 
site Coulomb repulsion. It has been shown [1,7] that the average of two electron operators 
on neighboring lattice site is proportional to the band's kinetic energy, it affects the kinetic 
energy, and it can be calculated in the simple way. The results have decreased the magnitude 
of interaction constants necessary to create the ferromagnetism. 

We have also shown [9] that the simple rectangular and parabolic density of state (DOS) 
cannot represent the itinerant magnetic elements. The rectangular DOS does not have a 
unique relationship between the Stoner field and magnetization, and the parabolic DOS 
does not allow for taking into account the inter-site interactions. The simplest possible 
DOS, which can be used, is the semi-elliptic DOS. 

In this paper we will use this new many body technique (mentioned above) in combi- 
nation with the simple semi-elliptic DOS, for calculation of the Curie temperature of basic 
ferromagnetic elements. 

2. The Model 

We follow closely our previous paper [7] and we employ the following model Hamiltonian 

H = - E t^j (cf^Cj^ + h.C.) - /iO E ^i - -^ E Uaflia + f/ E fliafli^a 

<ij>a i ia icr /-i \ 



+V E ninj + J E cZrCpCi^'Cja + J' Y. {ct]-ctiCjiCj^- + h.C. 

<ij> <ij>ij,a' <ij> 

where fio is the chemical potential, cf^{cicr) creates (annihilates) the electron with spin a 
on the ith lattice site, Ujo- = c^Qo- is the particle number operator for electrons with spin 
a on the ith lattice site, hi = hia + hi^a is the charge operator, U is the on-site Coulomb 
repulsion and F is the intra-atomic Hund's interaction. In the Hamiltonian (1) we have 
three explicit inter-site interactions; J -exchange interaction, J' -pair hopping interaction, 
V -density-density interaction. The spin dependent correlation hopping t^- depends on the 
occupation of sites i and j, and in the operator form it can be expressed as 

t1^ = t(l - hi_^){l - hj^„) + ti [hi_„{l - hj_^) + hj_^{l - hi_^)] + t2hi^ahj^a, (2) 

where t is the hopping amplitude for an electron of spin a when both sites i and j are empty. 
Parameters ti, ^2 are the hopping amplitudes for an electron of spin a when one or both 
of the sites i or j are occupied by an electron with opposite spin, respectively. Including the 
occupationally dependent hopping given by Eq. (2) into the Hamiltonian (1) we obtain the 
following result 



H = - E [to- ^t (^i-a + ^i-a) + 2tea;hi^ahj-a] Uf^Cj^ + h.C.) - /iQ E ^i+ , , 

<ij>(r ^ ^ i , (3) 

+potential energy 



where 



At = t-tu t,^ = ^^-h. (4) 

In this form it is quite visible that the kinetic interactions: At and t^^ are also the 
inter-site interactions. 

We assume that ti/t = S and t2/ti = 5*1. In general these parameters are different and 
they both fulfill the condition S* < 1 and 5*1 < 1 which is equivalent to t > ti > t2 (see Ref. 
[10]). For simplicity we will assume now that S = Si . With this assumption we have the 
relationship 

At = t-h=t{l-S), te. = i±^-ti = ^. (5) 

In the Hamiltonian (3) there are many inter-site interactions: At, tg^, J, J', V, for 
which we will use the modified Hartree-Fock (H-F) approximation. The main point of this 
approximation is to retain the following inter-site averages: Jo- = ( c^Cjo-), in addition to the 

usual on site averages: Ua = (cf^Cicr), which contribute to the Stoner field. 
For the average kinetic energy of +a electrons {K'^) one can write 



{K'^) = -t:ff ^ {ctc,^} = -ztlffh = -DJ., (6) 

where z is the number of the nearest neighbors, D^ = zf^rr is the half band-width of the 
+a electrons, and effective hopping integral f^rr = t ■b'^ , with b'^ defined by eq. (10) below. 
The kinetic energy can be calculated straightforward as 

D 

{K'^) = J f{E'^{e))E'^{e)p„{e)de = -DJ,, (7) 

-D 

with E'^{e) given by eq. (9) with omitted index k. Eq. (7) means that we can calculate 
our average product of two nearest neighbor operators: I„ ~ {K'^). In the case of zero 
temperature, T = K, and constant DOS: p^- i^) = const = -^ for —D < e < D, we obtain 
from eq. (7) that 

I„ = n„(l - n^) (8) 

This approximation will be used later on throughout this paper. 

The modified H-F approximation will lead to the modified dispersion relationship (see 

[1,7]) 

El = Sub'' - aAE^, (9) 

with the parameter 6°" describing the spin dependent change of the bandwidth given by 
.At . Jex / 2 ^0 .. \ J-V^ J+J' 



b. = l- 2— n_^ + 2^ [n\ - l\ - IJ_,) —h —^— (1°) 

and the effective Stoner shift (coming from all interaction constants standing in front of the 
single site operator rijo-) equal to 



TTl 

AE^ = {F-zJ)- + 2zl_^ {t,^n„ -At). (11) 



3. The Correlation Effects 

As in the previous paper [8], we employ the coherent potential approximation to the 
strong on-site Coulomb repulsion U and in the result we obtain in Fig. 1 the following 
schematic transformation of the initial paramagnetic DOS (see [11]). 

Fie. 1 



After applying the modified H-F approximation to the inter-site interactions we obtain 
the DOS de-formed by the inter-site correlation as shown in Fig. 2 (see 

Fig. 2 



4. The Numerical Results 

The magnetization in Bohr's magnetons is given by 

m = na - n^a, (12) 

with the electron occupation numbers for spin ±o" given by the following equation (see [8]) 



oo 



ds 

n±. = J p±. (e) ^^^(,_^^^^^^)/,^ , (13) 

— oo 

where the spin dependent DOS, p±a {£), is de-formed by both the on-site and the inter-site 
correlation, as in Figs 1 and 2. 

4.1 We consider initially the case of on-site Coulomb repulsion U = 

(i) There are only inter-atomic interactions V, J and J' 

We assume temporarily that the kinetic interactions are zero: At = tgx = 0. In this case 
the total Stoner field (from eq. (11)) is given by 



tot 



F + zJ. (14) 



To limit the number of free parameters we assume also that J' = J, and V = 0. 

The results for the minimum on-site Stoner field F = Ftot — zJ are shown in Fig. 3. They 
show a strong decrease of the critical on-site Stoner field under the influence of exchange 
inter-site interaction. 



Fig. 3 



Next, we computed from eqs (12) and (13) the temperature at which the value of mag- 
netization drops to zero in the presence of only the on-site Stoner field F and the inter-site 
exchange interaction J . The results are collected in Table I, for the values of electron oc- 
cupation and magnetization corresponding to ferromagnetic 3d elements. One can see that 
with the growing component of the on-site field (growing column number from 1 to 4) the 
Curie temperature keeps increasing well out of the range of experimental values (column 
no. 5). 

TABLE I 



(ii) the case of nonzero kinetic interactions At = t — ti and tgx = ^^ — ti 
We will assume now that V = J = J' = Q. In this case the dependence of the on-site 
critical field Fcr on electron occupation has already been investigated in [7] . The results have 
shown a strong decrease of the on-site critical field, with the decreasing hopping 'inhibiting' 
factor S. This decrease will also cause the decrease of Curie temperatures, as is shown in 
Table II for different values of S. 

TABLE II 



(iii) both kinetic interactions and inter-atomic interactions are nonzero 

The kinetic interactions are represented by the hopping 'inhibiting' factor S*, and the 
inter-atomic interactions by J = J' (V = 0). We still neglect the on-site strong Coulomb 
repulsion, f/ = 0. The results are shown in Fig. 4. It is evident from this figure that both 
S, and J strongly decrease the on-site Stoner field required to create ferromagnetism. 

Fig. 4 



Therefore we calculated the Curie temperature under the influence of these two factors. 
The results are shown in the Table III. 

TABLE III 



We can see from Table III that including simultaneously the kinetic interactions {S) and 
the inter-site exchange interactions (J) decreases the Curie temperature close to experimen- 
tal values (compare columns no.l and no. 5). 

4.2 The case of on-site Coulomb repulsion U ^ D 

The results in this case without the inter-site interactions are shown in Fig. 5. As 
can be seen from this figure there is very little difference in F^r for the realistic U = 3D, as 



compared to the case oi U ^ D. Therefore, in further calculations, we will use for simplicity 
large U :^ D. 

Fig. 5 



We will assume that both kinetic and inter-atomic interactions are nonzero. As previously 
these inter-site interactions will be represented by parameters S and J. The critical on-site 
Stoner field calculated in this case is shown in Fig. 6. The band is being split by the 
strong Coulomb repulsion, see Fig. 1 for the schematic depiction of the band. This is why 
the curves with S ^ 1 have a discontinuity at half filling, when the Fermi energy jumps 
from the lower to the upper sub-band. Comparing this figure with Fig. 4 shows that the 
influence of strong Coulomb on-site correlation is not really beneficial for n > 1.4, where 
the ferromagnetic 3d elements are located. 

Fig. 6 



To analyze further the influence of U on the magnetic properties we have calculated 
the magnetization versus temperature at f/ = and f/ = oo for the electron occupation 
representing iron, see Fig. 7. 

Fig. 7 



This figure shows that U does not decrease the Curie temperature at concentrations 
corresponding to existing ferromagnetic elements. Table IV shows the results of the Curie 
temperature in the presence of the same interactions as in Table III, with the difference of 
adding to it a very strong interaction of U. 

TABLE IV 



Comparing Table III with the Table IV tells us that including the on-site Coulomb 
correlation in the CPA approximation increases slightly the Curie temperature pushing it 
further away from the experimental values. 



5. Conclusions 

First, the direct calculations within the original Stoner model (see [13,14]) have pointed 
out that the Curie temperature, after adjusting the Stoner shift to the experimental magnetic 
moment at zero temperature, is much too high (see Table I, column 6 in here). 

As we have found previously [8] the inter-site interactions create the correlation effect (see 
Fig. 2) already in the first order modified H-F approximation. In this paper we introduced 
the kinetic and the inter-site interactions in the first order modified approximation that 

6 



allowed us to obtain the correct Curie temperatures and simultaneously the correct magnetic 
moments at zero temperature for 3d transition elements (see Table III). Considering the 
simplicity of the model, the results are close enough to the experimental data (see Table 
III). One can write symbolically that 



« 



Fij (m) Fi^i (m) ' 

where Fi^j{m) is the inter-atomic Stoner field fitted to the magnetic moment at zero tem- 
perature, and Fi^i{m) is the on-site Stoner field fitted to the magnetic moment at zero 
temperature. Intuitively speaking we can say that the inter-atomic field is much 'softer' 
and decreases faster with the temperature than the on-site Stoner field used in the original 
Stoner model (see Hubbard [15]). 

Hubbard's [15] calculations 'imply that two energy scales are operative in iron, one 
of the order of electron volts which is characteristic of the itinerant behavior (e.g., the 
bandwidth and the exchange fields), and another of the order of one tenth of an electron 
volt characteristic of the "localized" behavior (e.g., ksTc, the AE(y)).^ In our model this 
larger field would be the on-site field creating local moments existing even above the Curie 
temperature, and the smaller field would be the inter-site field responsible for their ordering. 

Adding up the on-site strong correlation U in the Section 4.2 (even in the CPA approxi- 
mation equivalent to the sum of some infinite set of diagrams (see [11])), has decreased the 
critical field, and enhanced magnetic susceptibility (see [16]), but did not improve the values 
of the theoretical Curie temperature (see Table IV). This fact can be understood better 
after examining Fig. 6, where the critical field (initializing magnetization) has dropped at 
half-filling but not at the end of the band, where the 3d elements are located. As was already 
established (e.g. see [17,18]) the electron correlation U can help in creating antiferromag- 
netism (AF) at the half-filled band, where the antiferromagnetic 3d elements (Cr, Mn) are 
located, by reducing the critical field for AF to zero. As was mentioned above, the model is 
very simple, the details of the realistic DOS could be included as well as the magnetization 
decrease through the spin waves excitation, which would bring the theoretical results to 
complete agreement with the experimental data. 
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Figure Captions 

Fig.l. Schematic DOS showing the influence of the strong on-site Coulomb correlation, 
U. The paramagnetic DOS for both spins, ±(T are solid lines. When U is strong enough, 
the band is split into two sub-bands, lower sub-bands have the capacity of 1 — n^^j for +a 
electrons, and l — ricr for —a electrons. The changes in the spin electron densities integrated 
over energy (to the Fermi level) are shown as the shaded areas in this figure, and they 
contribute with the sign to the correlation factor Ku. 

Fig. 2. Schematic DOS shows the influence of the inter-site interactions. The paramagnetic 
DOS for both spins, ±a are solid lines. The inter-site interactions change the relative width 
of the bands with respect to each other (described by the b'^ factors, see eq. (9)). The Stoner 
field, which would displace the bands with respect to each other, is assumed to be nonzero. 
The shaded areas in this figure are the inter-site correlation factor K^j. 

Fig. 3. Dependence of the on-site critical Stoner field, Fcr = F^^^ — zJ (in the units of half 
band- width), on electron occupation, for different exchange inter-site interactions J; J = 
(Stoner model) - continuous line, J = 0.5t - dashed hue, J = 0.75t - dotted hue. 

Fig. 4. Dependence of the on-site critical Stoner field, Fcr (in the units of half band-width), 
on the electron occupation, for different values of S and SJ; S* = 1 and J = (the Stoner 
model) - dot-dashed line, S = 0.6 and J = 0.5t - solid line, S* = 1 and J = 0.5t - dashed 
line, S = 0.6 and J = - dotted line. 

Fig. 5. Dependence of the on-site critical Stoner field, Fcr (in the units of half band- width) 
on electron occupation, for different on-site Coulomb repulsion U {J = 0, S = 1). 

Fig. 6. Dependence of the on-site critical Stoner field, Fcr (in the units of half band-width), 
on electron occupation, shows the influence of inter-site interactions, J = J', V = 0, and 
of the hopping interactions. At, tex, represented by the hopping 'inhibiting' factor S, in the 
presence of strong on-site correlation, U = oo. 

Fig. 7. Magnetization (in Bohr's magnetons) dependence on temperature, fn{T), with and 
without the strong Coulomb correlation, n = 1.4, D = 2.8eV for different values of parame- 
ters U, J and S; U = 0, J = 0, S = 1- solid line, U = oo, J = 0, S = 1 - dashed line, U = 0, 
J = 0.5t, S = 0.6 - dotted hue, U = oo, J = 0.5t, S = 0.6 - dot-dashed line. 
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TABLES 



El. n m D [eV] 



Tc[K] 



No. 1 


No. 2 


No. 3 


No. 4 


Ftot = 0, 


Ftot = zj{m) 


Ftot = zj + F{m) 


Ftot = 


F = -zJ(m) 


F = 


J = 0.75t 


J = 



No. 5 No. 6 



Fe 1.4 


0.44 


2.8 


1750 


3340 


3640 


4290 


1043 


4820 


Co 1.65 


0.344 


2.65 


1590 


3000 


3500 


4710 


1388 


3380 


Ni 1.87 


0.122 


2.35 


920 


1380 


1620 


1960 


627 


1020 



TABLE L Curie Temperature for Ferromagnetic Elements, the role of inter-site ex- 
change interaction, n=3d electron number normalized to 2, m=niagnetic moment at OK normal- 
ized to 2, D= 3d half band-width according to [12], columns no. 1,2,3,4- Tc[K]=Curie temperature 
for different are for different F and J, F{m) and J{m) are fitted to m at OK, T^'-^P experimental 
Curie temperature in K, T|'"^"^=Curie temperature for the Stoner model with parabolic DOS. 
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El. 


n 


m 


Tc[K] 








5 = 0.2 


5 = 0.6 


S = 1 


rr^ 


Fe 
Co 

Ni 


1.4 
1.65 

1.87 


0.44 

0.344 

0.122 


2445 
1580 
220 


3295 
3300 

870 


4290 
4710 
1960 


1043 
1388 
627 



TABLE II. Curie Temperatures for Ferromagnetic Elements, the role of kinetic 
interactions Tc = Tc{S) 



bl. n m 



No. 1 


No. 2 


No. 3 


No. 4 


No. 5 


5 = 0.6 


5 = 0.6 


5=1 


5=1 


T^P 


J = 0.5t 


J = 


J = 0.5t 


J = 





Fe 


1.4 


0.44 


2050 


3295 


3980 


4290 


1043 


Co 


1.65 


0.344 


1690 


3300 


3880 


4710 


1388 


Ni 


1.87 


0.122 


620 


870 


1720 


1960 


627 



TABLE III. Curie Temperatures for Ferromagnetic Elements, the role of kinetic and 
inter-site exchange interactions at [/ = 0, no. 1,2,3- Tc{S^ J)= Curie temperature for fixed J, S 
and F{m), F{m) is fitted to m at OK, nO. 4 -Tc{S = 1, J = 0)= Curie temperature for the case of 
Stoner model, no inter-site interactions {J = J' = V = At = tex = 0) , F{m) is fitted to m at OK 
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El. n m 



TclK] 



No. 1 


No. 2 


No. 3 


No. 4 


No. 5 


5 = 0.6 


5 = 0.6 


5=1 


5=1 


T^P 


J = 0.5t 


J = 


J = 0.5t 


J = 





Fe 


1.4 


0.44 


2420 


4000 


3810 


4960 


1043 


Co 


1.65 


0.344 


1850 


3270 


3990 


5180 


1388 


Ni 


1.87 


0.122 


630 


910 


1770 


2020 


627 



TABLE IV. Curie Temperatures for Ferromagnetic Elements, the role of kinetic and 
inter-site exchange interactions at t/ = oo 
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